Secular variability, geodetic precession and moment of inertia of binary pulsars 
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More and more binary pulsars show significant secular variations, in which the measured projected 
semi-major axis, x obs , and the first derivative of orbital period, Pb° bB , are several order of magnitude 
larger than the prediction of general relativity (GR). This paper shows that the geodetic precession 
induced orbital effects can explain both x and Pb measured in binary pulsars. Moreover, by this 
model we can automatically estimate the magnitude of the spin angular momenta of the pulsar and 
its companion star, and therefore the moment of inertia ( 10 44 g cm 2 to 10 45 g cm 2 ) of pulsar of binary 
pulsar systems, which agrees well with theoretical predictions. In other words, the contamination 
(residuals represented by x obs and P b obs ) in pulsar timing measurements might be caused by geodetic 
precession, an interesting gravitational effect we have been seeking for. 

PACS numbers: 04.80.-y, 04.80.Cc, 97.10.-q 



I. INTRODUCTION 

In the first gravitational two-body equation including 
spins 0, the precession velocity of the orbital angular mo- 
mentum vector, L, was expressed as around a vector, 
combined by the spin angular momenta of the pulsar, 
Si, and its companion S2. This velocity has been ig- 
nored, since it is much smaller, 2 Post-Newtonian order 
(PPN), than the precession velocity of the pulsar spin 
axis, 1.5PPN in binary pulsar systems. 

However, the small orbital precession velocity relative 
to the vector combined by Si and S2 doesn't mean that 
it is also small relative to the line of sight (LOS). Then in 
the study of observational effect of the orbital precession 
due to the coupling of the spin induced quadrupole mo- 
ment with the orbital angular momentum [2J. y. |j] . the 
precession of L is expressed as relative the total angular 
momentum vector, J, which is at rest to LOS (counting 
out the proper motion effect); rather than the vector (Si 
and S2), which are precessing relative to the LOS rapidly 
themselves. 

Thus, there seems a contradictory between Barker and 
O'Connell's equation and that of other authors mention 
above on which direction L should precess around. Actu- 
ally they can be consistent in the following scenario 0,0- 
In which L, Si and S2 all precess around J rapidly 
(1.5PPN), while the relative velocities of L relative to 
Si and S 2 are very small (2PPN). 

The coupling of the spin induced quadrupole moment 
and the orbit is 2 PPN p], so the quadrupole moment 
of the companion star has to be very large to make the 
effect observable. Therefore, this model is suitable for 
such binary pulsars, as neutron star-main sequence star 
(NS-MS), and neutron star-black hole (NS-BH) @ . 

Unlike NS-MS and NS-BH binaries, neutron star-white 
dwarf (NS-WD) and neutron star-neutron star (NS-NS) 
binaries have much smaller quadrupole moment, and the 
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spin angular momenta of the two stars can be compara- 
ble. To explain the observations of such binaries, alter- 
native model is necessary. 

Apostolatos et al[(J and KidderQ studied the preces- 
sion of L as relative to J in discussing the modulation of 
gravitational wave by the spin-orbit (S-L) coupling in co- 
alescing binary systems (actually NS-BH) . Similar to the 
former quadrupole-orbit coupling, this S-L coupling only 
considers the spin of the companion star with the orbit, 
and the spin of the pulsar is ignored. The difference is 
that the former quadrupole-orbit coupling corresponds to 
2 PPN, whereas the latter S-L coupling corresponds to 1.5 
PPN. Therefore, the latter is easier to explain the mea- 
surements of NS-WD and NS-NS binaries, which have 
smaller quadrupole moments and comparable spins. 

This paper develops the equation of orbital precession 
of Apostolatos et al and Kidder, making it applicable for 
general binary pulsars. We point out that once a binary 
system has two spins instead of one (the masses of the 
binary system can be different), the precession of the 
spin angular momenta of the two stars, Si and S2, lead 
to a variable S (S = Si + S2), which in turn tilts the 
orbital plane. Thus the additional motion of the orbital 
plane can explain both x and Pb measured in NS-WDs, 
i.e., PSR B1957+20 8] and PSR J2O51-O8270, which 
have been interpreted separately by different models 0, 
IToj ]. Moreover, the derivatives, Pb, Pb, and x can also be 
naturally interpreted. 

The geodetic precession induced x and Pb include an 
angle, X^j (the angle between L and J), which represents 
the intensity of S-L coupling (Xlj ~ S/L). Comparing 
with the measured x ohs , the magnitude of S/L can be ob- 
tained. Since L can be obtained from the binary parame- 
ters measured, we can thus estimate the magnitude of Si 
and S2, and in turn the moment of inertia of the pulsar, 
through the measured pulsar period. The obtained mo- 
ment of inertia is about 10 44 g cm 2 and 10 45 g cm 2 , which 
is in the range of strange stars and NSs. 

On the other hand, the geodetic precession model is 
supported not only by the secular variabilities measured, 
but also by theoretical predictions of the moment of in- 
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ertia of pulsars. 

Section II introduces the orbital precession in special 
cases. Section III discusses the orbital precession in gen- 
eral cases and its relationship with the orbital precession 
in special cases. Section IV applies the general orbital 
precession to two NS-WD binaries, explains the signifi- 
cant secular variabilities measured and estimates the mo- 
ment of inertia of pulsars in the two binaries. 

II. ORBITAL PRECESSION IN SPECIAL CASES 

The motion of a binary system can be seen as the pre- 
cession of three vectors, the spin angular momenta of the 
pulsar and its companion star, Si and S2, and the orbital 
angular momentum L. The change of the orbital period 
due to the gravitational radiation is 2.5 PPN; whereas 
the geodetic precession corresponds to 1.5 PPN. So the 
influence of gravitational radiation on the motion of a bi- 
nary system can be ignored in the discussion of dynamics 
of a binary pulsar system. Therefore, the total angular 
momentum, J = L + Si + S 2 , can be treated as invariable 
both in magnitude and direction (J = 0). With fio, the 
precession rate of L around J, the S-L coupling can be 
expressed asQ 

n x l = -fti x Si - n 2 x s 2 , (1) 

where fii and f2 2 represent the precessions of the pul- 
sar and its companion star, respectively. Ignoring those 
terms that are over 2 PPN, £!i and Vl 2 can be written 
asQ 

_ L 3m 2 . „ L 3mi. 

r! i = T^( 4 +— ) - ^2 = ^(4+— ) , (2) 
2r° mi Zr° m 2 

where mi and m 2 are masses of the pulsar and the com- 
panion star, respectively, and r is the separation of mi 
and m 2 . Notice L oc r 1 / 2 , Sli and f2 2 are 1.5 PPN. 

In the case of one spinning body, i.e., S 2 = 0, S = Si, 
L and S precess about the fixed vector J at the same rate 
with a precession frequency approximately 7J , 

I J I . , 3m 2 . ,„. 

Eq© is also correct if the two bodies have equal masses 
(Sli = fl 2 in EqlJTJ). 

Section III shows that the orbital precession in a gen- 
eral binary pulsar system (S2 ^ 0, S2 ^ 0, mi ^ m 2 ) 
will automatically lead to the significant variability in 
and x. 

III. ORBITAL PRECESSION AND ITS 
EFFECTS IN GENERAL CASES 

EqQ indicts that at any instance the variation of the 
angular momenta L (left) equals the variation of the an- 
gular momenta Si and S 2 (right). In magnitude, the 



right-hand side of EqQ is (recall S = Si + S 2 ) 

(ti) r = Cl 2 SsmX L s + (fii - £l 2 )Sl sinA L Si ■ ( 4 ) 

L can react to the torque, (ti)r, by precessing around 
J, with a very small opening angle of the precession cone 
Al.7 ~ S/L. Since the triangle (L, J and S) constraint 
much be satisfied at any instant. Thus the left-hand side 
of Eq|[[J can be written as 

(ti)l = |O nj X Ln | = CIoLsixiXlj , (5) 

where nj and no are unit vectors of J and L, respectively. 
By Eq(0} and EqlJSJ, we have the precession rate of L 
around J, 

fl = fl 2 sm\ LS + (Qi - n 2 )-g-sm\ LSl , (6) 

where sf — 5*1 cos 77551, which denotes the component 
of Si in the plane determined by S and J. Note that 
LsinAij « S is used in Eq@, since S/L <C 1. The 
difference between Eq© and Eq@ is obviously due to 
different number of spins included. Notice that the right- 
hand side of Eq© can as well be written by replacing 
subscribes 1 with 2 and 2 with 1. 

Thus the orbital precession velocity of Eq©, Slo, is de- 
rived via two assumptions: the conservation of the total 
angular momentum; and the triangle constraint on the 
precession cone of the orbit, Xlj- 

Notice that n Q of Eq© is 1.5 PPN. Which can be ab- 
sorbed by the advance of the precession of the periastron, 
to. The measured Co is givenUBl by: 

uj obs = uj gr + n cos Xlj « u GR + n . (7) 

Since Si and S2 precess with different velocities, fii and 
2 respectively (mi ^ TO2), S varies both in magnitude 
and direction (Si, S2 and S form a triangle). From the 
triangle of S, L and J, in reaction to the variation of S, 
L must vary in direction (|L| —const), which means the 
variation of Xlj (J is invariable). Therefore, by EqQ, a 
variable Lu obs is expected. 

The change of Xlj also means the orbital plane tilts 
back and forth, as shown in Fig 1, in turn both Xls and 
Xjs vary with time. Therefore, by Eq©, the derivative 
of the rate of orbital precession can be given by, 

n = n 2 Qi 2 x 3 X4 — f2i2Xi(fioi^2 + ^12^3) ) (8) 

where f2i 2 = Clx — fl 2 , f2oi = fii — Oo, X% — sin Xlsh 

X 2 = tann„i, X 3 = Svl ^ V2 sinT '\ l32 , and X A = cos \ Xls ; 

with a = sin A./5 + c °^^ , S V1 = S 1 sin A JS1 and S V2 = 
S' 2 sinAj5 2 represent components of Si and S 2 that are 
vertical to J. 

Note that fli and £l 2 are unchanged, and Xls are 
unchanged, since they decay much slower than that of 
the orbit |g] ■ The second derivative of Slo are given by, 

f2o = ^2^12(^3^4 + X3X4) — Qi 2 Xi (f2oi AT 2 + Qi 2 X 3 ) 
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-fii2Xi(fioiX 2 + niaX 3 ) , (9) 

where X a is the first derivative of X a . £7 can be easily 
obtained from Eq@. 

Ooj the derivative of FIq, can be absorbed by Pb. The 
variation in the precession velocity of the orbit results in 
a variation of orbital frequency (i/j = 2n/Pb), v' h — Vb = 
O A.t, then we have i>b — f^o an d Vb = ^Oj therefore, 

b ~ 27T ' b ~ 2tT 7T ' [ ' 

From Eq(|SJ and Ea illOfl . we can see that the contribution 
of the orbital precession to Pb is 1 PPN, which is much 
larger than the contribution of GR to Pb- Since it de- 
pends on parameters, Xi, X 2 , X 3 , and Xls of Eq©, then 
it is also possible that the geodetic precession induced Pb 
is comparable to, or even smaller than the prediction of 
GR in special cases (special combination of Als, Xlsi, 
■qssi and Si/S). 

The effect of S-L coupling on secular evolution of the 
orbital inclination, i, can be given as, 

cos i = cos Xlj cos I — sin Xlj sin /cos 770 j (H) 

where / is the angle between the total angular momen- 
tum, J, and LOS, and rjQ = f2 £ + r\i is the initial 
phase) is the phase of precession of L. Thus i is also a 
function of time. By Ea (|ll|l we have, 

di/dt= -^sin?7o + 0(r- 7 / 2 ) . (12) 

Thus all parameters that are related with the orbital in- 
clination will change with the geodetic precession. Hence 
the projected semi- major axis, x = aisinz/c, will vary 
as the precession of the orbital plane. By Ea (|12|l . the 
derivatives of the projected semi-major axis are 

x = — xflo sin Xlj sin^o cot i , (13) 



IV. EFFECTS OF GEODETIC PRECESSION 
INDUCED Q 

The S-L coupling induced flo leads to the precession of 
the orbital plane, then the orbital inclination, z, varies, 
as in Ea l|12|l . and in turn the projected semi-major axis, 
x, varies, as in Ea i|13fl . Comparing the predicted x 
with the observational x obs , we can constrain the angle 
Xlj ~ S/L, and therefore the moment of the inertia of 
the pulsar. 



A. PSR J2051 0827 

The measured orbital period derivatives of 
PSR J2O51-O8270 are list in Table I, the deriva- 
tives of the semi-major axis are x obs = —23(3) x 10~ 14 
and x obs = 7(15) x 10~ 22 s — 1 . The measured ones are 
much larger than the corresponding predictions of GR, 
pGR = _3.g x 1(T 14 and x GR = -1.3 x 1CT 19 , 

as well as the proper motion induced effects, 
±pm = 4 (2) x 10 -i7 n\ 

With mi = 1.40M©, m 2 = O.O3M and P b = 
0.09911d|3, we have the semi-major axis, a — 
(GM/i/j 2 ) 1 ' 3 = 7.08 x 10 10 cm, and orbital angular mo- 
mentum, L = v b ^a 2 (l — e 2 ) 1 / 2 = 2.17 x 10 50 g cm 2 s _1 , 
with fi the reduced mass. From Eq(J2Jl, we have fl 2 = 
3.25xl0- 9 s _1 , and from Eq©, |fi | ~ ^2 « 3xl0 -9 s -1 . 
With x = x obs = 0.045052, |i ofcs | = 2.3 x 10~ 13 , EqO 
becomes 

0.045x(3xl0^ 9 )sinA L ./|sm77oCoti| = 2.3xl0~ 13 . (15) 

Thus sin Al j j sin 770 cot i I es 2 x 10~ 3 . If | sin 770 cot -i | ss 

1 x 10 _1 , then Xlj ~ 2x 10 -2 , which means S/L ~ 

2 x 10 -2 . With L obtained above, we have Si w S « 
4 x 10 48 g cm 2 s _1 . Having the measured pulsar period, 
P — 4.50864ms, the moment of inertia of the pulsar is 
Ji w 3 x 10 45 g cm 2 . 



B. PSR B1957+20 



'x = — i(f2o/^o + Qq cot 770) . (14) 

As discussed between Eq{7J) and Eq©, the change of S 
leads to the change of fioj an d m turn Pb and Pb. If 
there is only one spin, then S is a constant, thus there 
will be only a static addition (Qq = const) to the apsi- 
dal motion, as shown in EqQ. In this case, the orbital 
precession is also static, which will only cause a variation 
in x, as shown in Ea l|13fl . whereas, Pb and Pb are not 
influenced at all (CIq = and Oq = 0). This is why the 
quadrupole-orbit model (also the one spin S-L coupling 
model) cannot explain P b and P b in PSR J2051-0827 
and PSR B1957+20. 



The negative orbital period derivative changes to posi- 
tive during the observation of PSR B1957+20@. The 
measured derivatives of the orbital period are shown 
in Table 2. The measured upper-limit to x is x obs < 
3 x 10 -14 . The measured ones are also much larger than 
the corresponding GR predictions, P b GR = —3.4 x 10~ 15 
and x GR = —6.1 x 10~ 21 , as well as the proper motion 
induced effect, x pm = 2.4 x lO" 16 ^. 

Similarly, with mi = 1.40M Q , m 2 = 0.025M© and 
P b = 33001. 91484s|g, we have a = 1.74 x 10 11 cm, L = 
2.82 x 10 50 g cm 2 s- 1 , Q 2 = 3.43 x lO -1 ^ -1 , and in turn 
|O | « tt 2 « 3xl0- 10 s -1 . FromEqlO, with x = x obs = 
0.0892253 and \x obs \ < 3 x 10" 14 , we have 

0.089x(3xl0" 10 )sinA L j|sin77 cot7:| < 3xl0~ 14 , (16) 
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which leads to , sin sin?7o cot i\ < 1 x 10 -3 . If 
| sin 770 cot i\ « 1 x 10 _1 , then X LJ < 1 x 10~ 2 . Which 
means S/L < 1 x 10~ 2 , and by the obtained L, we have 
Si ~ S < 3 x 10 48 g cm 2 s _1 . Finally, with the measured 
pulsar period, P — 1.607ms, the moment of inertia of the 
pulsar is h < 8 x 10 44 g cm 2 . 

The simple estimation above indicates that 
PSR B1957+20 has similar moment of inertia as 
that of PSR J2051-0827. 

In these two estimations, the assumptions |fio| ~ ^2 
and I sin 770 cot z| « 1 x 10 _1 are used. Consider the de- 
viation of these two assumptions from the corresponding 
true values, the most likely moment of inertia from the 
observational constraints is about 10 44 g cm 2 -10 45 g cm 2 . 

V. EFFECTS OF GEODETIC PRECESSION 

induced n , n , fio 

From Eq© and EqQ, we can see that tl , Cl and fio 
can vary in wider and wider range, which can be written 
approximately as (for the convenience of estimation), 

|o | « cin\ , |ii„| « C2 r» 3 , |n | » c 3 n 4 2 , (17) 

where ci represents the contribution of X a to Qo, C2 cor- 
responds to the contribution of X a and X a to f2oj an d 
C3 corresponds to the contribution of X a , X a and X a to 
tla- Obviously tlo can vary in a wider range than flo, 
and f2o can vary in a wider range than Oq • 

For PSR J2051-0827, ft 2 « 9.8 x 10~ 18 s- 2 , ft 3 , w 
3.2 x 10~ 26 s~ 3 and also (equation not shown) w 
1.0 x 10~ 34 s~ 4 . Assuming ci = c 2 = c 3 = 0.2, and 
by Ea l(T7)l and EaltTUIl. the derivatives of P b can be ob- 
tained, \P b \ k 2 x 10- 1X , \P b \ « 8 x lO-^s" 1 and 
|Pb| « 2 x 10 _28 s -2 . Which can be well consistent with 
the corresponding observations, as shown in Table I. 



For PSR B1957+20, the measured derivatives of P b obs is 
larger than that of PSR J2051— 0827. By assuming ci = 
1, C2 = 1 x 10 3 and C3 = 1 x 10 4 , and through the same 
treatment, we have \P b \ w 2 x 10 -11 , |P 6 | « 7 x 10~ 18 s _1 
and \ 'p b \ k 2 x 10~ 26 s -2 , which can be consistent with 
the corresponding observations, as shown in Table II. 

It is likely that in PSR B1957+20 the magnitude of Si 
and S2 are close, so that S can be small, which leads to 
a large X a , and in turn large derivatives of P b as shown 
in Eq© and Eq©-Eq|TO|). 



VI. DISCUSSION AND CONCLUSION 

Geodetic precession induced pulsar spin in binary pul- 
sar PSR B1913+16 has been studied through structure 
parameters of pulsar profile [I2T IT^ . Actually the varia- 
tion of pulsar spin axis (former) and the measured sec- 
ular variabilities discussed in this paper are induced by 
the same physics underlying in binary pulsar systems, S- 
L coupling. The former corresponds to the reaction of 
the pulsar spin to the torque in the S-L coupling; while 
the latter corresponds to reaction of the orbit to the same 
torque. 

By adding one more spin into the S-L coupling of Apos- 
tolatos et al and Kidder we establish the relation- 

ship among the secular variability, geodetic precession 
and moment of inertia of binary pulsars. The new model 
provides: (a) a unified model that explain both x and P b , 
which has been interpreted separately by different mod- 
els, (b) a new method which can extract the moment of 
inertia through pulsar timing measurement, (c) a new 
test of the geodetic precession which has strong effects to 
pulsar timing. 
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TABLE I: Measured parameters compare with the geodetic 
precession induced ones in PSR J2051— 0827 

observations predictions 

± obs ~ „23(3) x lCT 14 \x\^2x 10~ 13 a 

(x/±) ob3 < -3.0 x 10 _9 s- 1 \x/x\ » 3 x HTV 1 b 

P£ bs = -15.5(8) x 10" 12 |A| = |%^H ~ 2 x KT 11 

P£ bs = 2.1(3) x lO^V 1 \P b \ » |^| Ri8x lO^s" 1 

Pg bs = 3.6(6) x 10- 28 s~ 2 \P b \ « g 2 x 10- 28 s- 2 



We assume ci = C2 = C3 = 0.2 in Ea ll7H . a \x\ is given by 
Ea ll5l . and b x/x is given by Eg 1141 . 



TABLE II: Measured parameters compare with the geodetic 
precession induced ones in PSR B 1957+20 

observations predictions 

\x\° ba < 3 x HT 14 |i| < 3 x 10" 14 a 

Pg u = 1.47(8) x lO" 11 \P b \ = |^|^| Ri2x lO" 11 

P£ bs = 1.43(8) x HT 1 ^- 1 |A| ~ |%^| ~ 7 x lO" 1 ^" 1 

\P b \ obs < 3 x !Q- 26 s- 2 I Pj, I g |%^| w 2 x !Q- 26 s- 2 



We assume ci = 1, C2 = 1 X 10 3 , and C3 = 1 x 10 4 in Eg ATI. 
a \x\ is given by E0 II6II . 



FIG. 1: The orbital plane, n, tilts back and forth ( correspond- 
ing to L moves along the curves). L, S and J form a triangle, 
in which J is invariable both in magnitude and direction; L 
varies in direction only; and S varies both in magnitude and 
direction 



